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Quantum spin Hall insulator/metal interfaces are naturally formed in graphene ribbons with 
intrinsic spin-orbit coupling by selectively doping two regions creating a potential step. For a clean 
graphene ribbon, the transmission of the topological edge states through a n-n or p-p junction is 
perfect irrespective of the ribbon termination, width, and potential step parameters due to the 
orthogonality of incoming and outgoing edge channels. This is shown numerically for an arbitrary 
crystallographic orientation of the ribbon and proven analytically for zigzag and metallic armchair 
boundary conditions. In disordered ribbons, the orthogonality between left- and right-movers is in 
general destroyed and backscattering sets in. However, perfect transmission is restored by increasing 
the ribbon's width, even in the presence of strong edge roughness. 

PACS numbers: 73.43.-f, 73.23.Ad, 72.25.Mk, 73.40.-c 



I. INTRODUCTION 

Since the discovery of the quantum Hall effect^, edge 
states have played an important role in condensed mat- 
ter physics. With the prediction of the topological insu- 
lator phase^ ^ and its observation in HgTe/CdTe quan- 
tum wells^ ^, the study of edge state physics has received 
renewed attention, both from the theoretical and experi- 
mental side. This new state of matter is identified by an 
energy gap in the bulk with gapless states propagating 
around the edges of the sample. In contrast to the nor- 
mal quantum Hall effect, these edge states preserve time- 
reversal symmetry and must be spin-polarized as a conse- 
quence: at the system's boundary a pair of counterprop- 
agating spin-up and spin-down states exists as sketched 
in Fig. [ij As equal-spin states with opposite velocity are 
thus spatially separated at opposite edges of the sam- 
ple, and since different-spin states at the same boundary 
form a Kramers' doublet!^, edge states are protected 
against backscattering from time-reversal invariant im- 
purities. This so-called quantum spin Hall (QSH) state 
has been first studied theoretically in graphene^, where 
it results from the intrinsic spin-orbit coupling induced 
by the carbon atoms. 

Graphene inherits interesting properties from its linear 
dispersion for low-energy quasi-particle excitationJ^, 
mak ing i t a candidate for future spintronics applica- 
tion^SI^. Part of its functionality comes from the ability 
to tune the carrier density and carrier type (electron or 
hole) with electric gates, which allows for the straightfor- 
ward creation of n-p junctions in graphene. Such junc- 
tions allow the investigation of the dynamics of the effec- 
tive massless Dirac-fermions in graphene, showing inter- 
esting properties such as Klein tunneling^^ or beam 
collimation and lensing^^. Graphene's functionality can 
be extended by forming junctions with different mate- 
rials, such as ferromagnetJ^^l^ or superconductorJ^^'^. 
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FIG. 1: (Color online) a) A graphene ribbon with intrin- 
sic spin-orbit coupling is a topological insulator with gapless 
spin-polarized edge states at low energies. In the presence 
of a potential step of height Uo bigger than the spin-orbit 
gap Aso, a quantum spin Hall insulator/metal interface is 
formed where topological edge channels evolve into multiple 
bulk propagating modes, b) Energy bands for a zigzag termi- 
nated ribbon at the left and right sides of a potential step of 
height Uo and width 6. The energy range of interest is marked 
in grey. 



With the advent of topological insulators, the investiga- 
tion of n-p junctions in the QSH phase (in graphene 
and in quantum wellJ^ as well as interfaces between this 
phase and other material^^^M^ has just been initiated. 

In this paper we study the transmission properties of 
the QSH edge states in graphene ribbons through an 
electrostatic potential step, defined for example by se- 
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lectively doping left and right regions of the ribbon by 
means of two back gates. When the potential step height 
Uq is such that the Fermi energy in both regions lies in- 
side the spin-orbit induced band gap Aso, perfect trans- 
mission is trivially expected since transport in both re- 
gions occurs in the form of topological edge channels. 
Spin-preserving backscattering processes then amount to 
crossing the sample width and are forbidden as long as 
the ribbon's width is larger than the width of the edge 
states^^. However, when Uq is bigger than Aso, a QSH 
insulator/metal junction is formed and topological edge 
channels in one region evolve into multiple bulk channels 
above the bandgap in the other (see Fig. [T]). Backscat- 
tering of the QSH edge states could then in principle 
occur through the metallic region where a path opens 
up allowing the carriers to cross the sample's width from 
one boundary to the other. This reflection process would 
lead to a spin-current parallel to the interface (a similar 
spin-current has been predicted at an interface between 
semi-infinite graphene planes in the presence of Rashba 
spin-orbit coupling^^). However, we show in this paper 
that the topological edge states in a clean ribbon of an ar- 
bitrary crystallographic orientation cannot be reflected at 
a QSH insulator /metal junction (either n-n or p-p) due to 
the orthogonality of left- and right-moving edge states at 
opposite boundaries. A similar mechanism was recently 
shown to be responsible for the perfect transmission of 
the zero Landau level edge state at a magnetic flux step in 
folded graphene nanoribbons^^. Perfect transmission is 
also present in n-p junctions, but only when the ribbon is 
in the armchair or antizigzag configuration, because only 
then does the n-p junction not introduce intervalley scat- 
tering^. The orthogonality between edge states is proven 
by deriving analytic expressions for the QSH edge wave- 
functions in zigzag and armchair terminated ribbons. For 
general ribbon terminations, the absence of backscatter- 
ing is furthermore validated by numerical calculations. 
Interestingly, we also find that the topological edge state 
width depends on the ribbon's termination. For an arm- 
chair ribbon, this width only depends on the magnitude 
of the spin-orbit induced gap, whereas it is practically 
independent of this parameter for zigzag ribbons where 
the edge state width becomes energy-dependent instead. 
This observation is further confirmed by numerical sim- 
ulations of the charge density across the sample. These 
results for clean junctions are analyzed in Section [ll| and 
in the Appendix. 

The influence of disorder, produced for example by 
edge roughness, is analyzed in Section [TTTj In general, the 
orthogonality between equal-spin left- and right-movers 
is then destroyed due to intervalley scattering and the 
conductance of the junction decreases as inter-boundary 
backscattering becomes possible through the metallic re- 
gion. Backscattering is also present for intrinsic sources 
of intervalley scattering, like a (clean) n-p junction in a 
zigzag ribbon^^. However, we find that in any case per- 
fect transmission is restored upon increasing the ribbon's 
width ly, essentially because the probability of forward 



scattering increases with the ribbon's width while the 
only available backscattering channel remains one. The 
conductance recovery with increasing W depends on the 
ribbon termination due to the different edge state char- 
acter for different ribbon boundaries. A summary of all 
our results and some concluding remarks are given in 
Section IVl 



II. EDGE STATES IN A BALLISTIC JUNCTION 

A. System description 

We consider the system depicted schematically in 
Fig.[l] It consists of a graphene ribbon subject to the in- 
trinsic spin-orbit coupling which opens a gap in the bulk 
graphene spectrum. The Fermi energy in the left part of 
the ribbon is considered to lay inside this gap, so that 
transport in this region can only take place through the 
topological edge channels (see Fig. [T]). In the right part 
of the ribbon, the potential is shifted upwards or down- 
wards by an electrostatic gate in order to reach doping 
levels above the spin-orbit induced bulk gap. Topologi- 
cal edge states then do not play a role anymore in this 
region and transport is mediated by normal propagating 
bulk modes. In this way we create a QSH insulator/metal 
junction in graphene (although the spin-orbit interaction 
is present throughout all the ribbon). The Hamiltonian 
of our system is given by 

H = Ho^Hso^H^. (1) 

The first term represents the standard hopping between 
nearest-neighbor carbon atoms, 

Ho = -t ^ 4,^9,,^, (2) 

where ^ is the creation operator for an electron with 
spin Sz on site i. The second term in the Hamiltonian 
accounts for the intrinsic spin-orbit coupling introduced 
by Kane and Mele^: 

Hso = itso Yl 4,5. ^j^s, ' (3) 

where the summation runs over next-nearest neighbors 
(note that in Eq. ([2| hopping occurs between different 
sublattices A and while in (|3| hopping is within the 
same sublattice). The prefactor = +1(— 1) if the elec- 
tron takes a left (right) turn to reach its next-nearest 
neighbor. Since this term commutes with the electron's 
spin the Hamiltonian of the system decouples into 
two independent Hamiltonians for spin-up and spin-down 
that can be considered independently. In this case, we 
have edge states for each boundary with spin-up moving 
in one direction and spin-down moving in the opposite 
one (see Fig. fl]). In the presence of spin-mixing terms 
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that do not break time-reversal symmetry (like for ex- 
ample the Rashba term) , conserved currents do not have 
a well-defined spin projection. Nevertheless, it will still 
be possible to define two sets of counter-propagating cur- 
rents at each boundary that are topologically protected 
against non-magnetic impuritie^^. For simplicity, below 
we consider the case where is a good quantum num- 
ber, but the physics we describe should be analogous if 
time-reversal symmetric spin-mixing terms were present. 

The last term H\j in the Hamiltonian describes the 
on-site graphene doping, given by 



(4) 



The sketch in Fig. [T] shows a n-n junction. In this case 
the doping profile along the boundary direction can be 
modeled by a smooth potential step of height Uq of the 
form 



U{x) = 



Uo 



1 +erf 



2{x-L/2) 



(5) 



where L is the length of the ribbon and S controls the 
smoothness of the step. This gives U{x) = when 
X — L/2 < 5 (left part of the sample), while U{x) = —Uq 
for X — L/2 > S (right part of the sample). In order to 
have a QSH insulator/metal junction, we need Uq > Aso 
and consider the energy range [0, Aso] depicted in grey 
in Fig. [ijb). Here Aso is the spin-orbit induced bulk 
gap (for an infinitely wide ribbon it is Aso = ^V^tsOi 
and slightly larger for finite width ones). An equiva- 
lent p-p junction is obtained by changing Uq by —Uq in 
Eq. (5|, and considering energies between —Aso and 0. 
Furthermore, we could consider n-p (or p-n) junction if 
we choose Uq > 2 Aso- In the following we will focus on 
the n-n junction shown in Fig. [T] for simplicity. 



B. Perfect transmission 

At a QSH insulator/metal junction in a graphene rib- 
bon, an incoming edge state may in principle either be 
transmitted or refiected at the junction. Spin-preserving 
refiection is only possible when the electron is transferred 
from the incoming edge state into the counterpropagat- 
ing edge state at the opposite boundary. Such an inter- 
boundary refiection is in principle possible because the 
electron's wavefunction to the right of the junction is not 
constrained to the edges, but spans the whole bulk in 
multiple propagating metallic channels. However, one 
of the main points of our paper is that, at a valley- 
preserving junction, this refiection process is forbidden 
due to the orthogonality of forward and backward topo- 
logical edge channels, much like in the Klein paradox. 
In Fig. [2] the conductance of the junction in units of 
Go = 2p/h is plotted as a function of the Fermi en- 
ergy of the incoming electrons (in units of t, the nearest- 
neighbor hopping amplitude). Different ribbon termina- 
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FIG. 2: (Color online) Conductance vs Fermi energy for a 
QSH insulator/metal n-n junction in the ballistic regime (see 
Fig. [1]). Perfect transmission is found irrespective of the 
ribbon termination: zigzag (solid red line), armchair (dashed 
blue line), and an intermediate crystallographic orientation 
(thin green line). All ribbons have a width W = 54a. The 
potential step height is Uo = 0.25t and its width S = 5a. 



tions are considered: armchair (parametrized by an an- 
gle 6> = 0, dashed blue curve), zigzag {6 = 7r/6, solid 
red curve), and an intermediate crystallographic orien- 
tation {0 = 0.47r/6, thin green curve). The ribbon's 
width is W = 54a and the potential step has param- 
eters U = 0.25t and 6 = 5a. Perfect transmission is 
manifested by a fiat first conductance plateau at energies 
where only edge states are present on the QSH insulator 
side. For higher energies, transport is carried by addi- 
tional propagating modes above the spin-orbit gap that 
are not protected by orthogonality and thus the trans- 
mission is reduced (compared to the situation without 
a potential step, dotted curves). Perfect transmission 
below the gap holds irrespective of the ribbon's crystal- 
lographic orientation and its width W,^^ and is indepen- 
dent of Uq > and 5 (as long as 5 is larger than a few 
carbon-carbon lattice constants). 

Our numerical results for mesoscopic transport are ob- 
tained using the recursive Green's function techniqu^^. 
A double-sweep generalization^^ gives us access to the 
local charge density in real space. This allows for a di- 
rect visualization of the shape of the edge states and 
their propagation across the potential step, as depicted in 
Fig. [3] for different ribbon terminations and energies be- 
low the bulk spin-orbit gap. The charge density is shown 
in arbitrary units, with black representing the absence of 
charge, and white its maximum value. Since spin-up and 
spin-down are decoupled in our problem, they can be an- 
alyzed independently and their contribution summed up 
at the end. In Fig. |3] only the spin current incoming from 
the upper boundary is considered. A similar contribu- 
tion (mirror symmetric with respect to the ribbon axes) 
coming in from the lower boundary with the opposite 
spin is not shown. In all ribbons, the incoming topo- 
logical edge state in the left region propagates towards 
the junction and then spreads out into multiple channels 
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FIG. 3: (Color online) Numerical results for the charge den- 
sity across a potential step {Uo = 0.25t and S = 5a) for a clean 
ribbon W = 54a with tso = 0.03t. A section L = 100a 
of the infinitely long ribbon is shown, and various ribbon ter- 
minations and Fermi energies below the spin-orbit gap are 
considered. At the left of the n-n junction (placed in the 
middle of each flake), transport is carried by a topological 
edge channel that evolves into multiple normal channels at 
the right. A similar contribution with opposite spin comes in 
from the lower boundary (not shown). 



available in the bulk of the metallic region. There is no 
sign of a backscattered edge state at the lower bound- 
ary on the QSH insulator side. It is interesting to note 
that the width of the topological edge state clearly de- 
pends on the crystallographic orientation of the ribbon: 
it is energy-independent for an armchair ribbon while it 
strongly depends on energy for the zigzag and interme- 
diate one. This behavior will have consequences in the 



presence of disorder, as analyzed in Section III 



For zigzag and armchair ribbons we were able to derive 
analytic expressions for the topological edge state wave- 
functions. For clean metallic ribbons, this allows us to 
prove the orthogonality of forward and backward chan- 
nels at opposite boundaries that is at the basis of the 
perfect transmission observed in Fig. [2] This orthogo- 
nality is in the valley sector of the wavefunction, which 
we assume must be preserved in the presence of a smooth 
electrostatic potential step. Furthermore, our analytical 
results also offer an explanation for the different edge 
state widths for armchair and zigzag ribbons observed 
in Fig. |3] Full details of the derivation are given in the 
Appendix. Below we present a summary of the most rel- 
evant results. 



C. Edge states at an armchair boundary 

First we consider an armchair ribbon extending along 
the vertical ^/-direction, so that the boundaries are at 
X = and x = W (see Fig. [9] in the Appendix). For 
simplicity, we consider that the ribbon is wide enough so 
that the two opposite boundaries can be treated indepen- 



dently. In the armchair case, it is known that the topolog- 
ical edge states have wavefunction solutions around the 
K and K' points in the first Brillouin zon^^. Expand- 
ing the Hamiltonian for small momentum q = {qx^Qy) 
around these points and considering small energies be- 
low the spin-orbit induced gap, we make the following 
ansatz for the edge states of a wide ribbon at the left 
boundary {x = 0): 



e '^^ aoe^^-^ 



This is a superposition of modes around the two inequiv- 
alent K and K' points {K'^ ^ = ±-^^^x) with amplitudes 
ao and Pq. These are propagating in the ^/-direction and 
exponentially decaying in the x-direction. |?/^) is a 

two-component spinor for the K (K') valley which con- 
tains the relative amplitude of the wavefunction on the A 
and B sublattices. In the Appendix we show that Eq. (|6| 
satisfies the armchair boundary condition when 



sgn((7^)sgn(£:). 



(7) 



where Sz = +1( — 1) for spin- up (spin-down) and sgn(g^), 
sgn(£^) are the signs of the particle's wave vector and 
energy, while Aso = "^V^tso is the spin-orbit induced 
gap for a wide ribbon. The dispersion of these states 
\^ E = ±.hvYqy^ where vp = 3at/2 is the Fermi ve- 
locity. As n needs to be positive for a bounded solu- 
tion, positive energy states (electron-like) either propa- 
gate with positive wavevector [sgn(g'^) = +1] and spin- 
down {sz = —1), or negative wavevector [sgn{qy) = — 1] 
and spin-up {sz = +1) along the left boundary. For the 
right boundary we use the same ansatz ([6| but now the 
ribbon extends to values x < W and the boundary con- 
dition is met ai X = W. Now needs to be negative 
and the positive energy edge state has positive (nega- 
tive) wavevector with spin-up (down). If we rotate our 
lattice by 90 degrees, we get the electronic edge states 
propagating along the sample's boundaries as sketched 
in Fig. [T] 

The width of the edge states is given by the parameter 
A = and is thus determined solely by the spin-orbit 

coupling strength. In our numerical calculations we have 
chosen tso = 0.03t (the same value used by Kane and 
Mele in Ref. 2, although this value is thought to be 
overestimated^^ ^^), giving a decay length A ^ 10a. This 
is in good agreement with the edge state width observed 
in Fig. [3| 

The orthogonality between forward- and backscattered 
edge states can be understood as follows. In the Ap- 
pendix, we show that the Hamiltonian of the armchair 
ribbon with spin-orbit coupling can be transformed into 
a valley-isotropic form,^ meaning that it has the same 
shape for both the K and valley. The valley degree 
of freedom (also called isospin) is thus conserved by the 
Hamiltonian. The potential step conserves the isospin as 
long as it is smooth on the scale of the lattice constant. 



5 



20 



15 



^ 10 



-armchair t, 
-zigzag t, 




FIG. 4: (Color Online) Comparison of the topological edge 
state width of armchair (closed black circles) and zigzag rib- 
bons (solid red line for A2 and open blue circles for Ai) 
as a function of Fermi energy. The spin-orbit coupling is 
tso — 0.03t. The dominant A2 contribution in the zigzag 
case is also plotted for different values of tso (0.01,0.03,0.06 
for dashed, solid and dotted red lines respectively). 



The isospin is described by a spinor |r) = (a,/3), where 
a and j3 are the coefficie nts of the plane waves in the K 
and K' valley [see Eq. ( |A20 )]. In the valley-isotropic 
basis, they are given by |r)^ = ^(is^,— 1) for the left 
boundary, and |r)^ = ^(is^e"^^^^, 1) for the right one. 



Writing W 



{N + l)a in terms of the number N 



of carbon atoms along the ribbon's width,SSl one can see 
that the isospins of counterpropagating equal-spin edge 
states (moving at opposite boundaries) are orthogonal 
when the armchair ribbon is metallic, i.e., for = 3n — 1, 
with n = 1, 2, .... Scattering between these states is thus 
forbidden which leads to the absence of reflection. 

The case of semiconducting armchair ribbons is a bit 
more involved. Let us show this with an example. Imag- 
ine we had a ribbon without potential step but with a 
smooth intrinsic spin-orbit step instead. The left region 
is then exactly the same as in Fig. [l] but the right region 
now is plain graphene (i.e., without spin-orbit coupling). 
For metallic armchair ribbons (actually for zigzag and in- 
termediate ones too), such a junction constitutes a QSH 
insulator/metal interface and we observe numerically a 
perfect transmission for energies below the induced spin- 
orbit bulk gap. However, the semiconducting armchair 
ribbon is insulating for small energies in the right region, 
and the transmission is zero for such energies. Therefore, 
knowledge of the wavefunction in the left region only is 
not sufficient in order to make a statement for the trans- 
mission in semiconducting armchair ribbons. In our case 
of a constant spin-orbit coupling and a potential step we 
can only state that numerically we also observe perfect 
transmission as long as the ribbon is metallic in the pres- 
ence of spin-orbit coupling, i.e., W > 2 A, even though 
the isospins of the incoming and outgoing edge states are 
not orthogonal. 



D. Edge states at a zigzag boundary 

Now we consider a wide zigzag ribbon that extends 
along the horizontal x-direction, so that the boundaries 
are ai y = and y = W (see Fig.[9]in the Appendix). For 
zigzag ribbons, the dispersive topological edge states are 
centered around the M point {qxCt = tt /V3). Therefore, 
the standard low energy expansion of the Hamiltonian 
around the K and K' points (that leads to the Dirac 
equation in the absence of spin-orbit interaction) is not 
reasonable anymore, since the low-energy solutions in the 
presence of the intrinsic spin-orbit coupling occur for mo- 
menta far away from those points. In Ref. 42, Zarea and 
Sandler resorted to a tight-binding method to calculate 
the wavefunctions and dispersion of the topological edge 
states. In the Appendix, we follow a slightly different 
approach that allows us to study the edge state width in 
some more detail. We find that the wavefunction at the 
lower boundary ^ = can be written as 



ae 



(8) 



It is interesting to note that one needs a superposi- 
tion of two plane waves in order to satisfy the boundary 
conditions, both for the armchair and the zigzag ribbon. 
In the absence of spin-orbit coupling, the wavefunction 
has to vanish at both A- and B- sites at the boundary 
of an armchair ribbon, but only at one sublattice for the 
zigag ribbon^^. However, as spin-orbit coupling induces 
next-nearest neighbor hoppings, the wavefunction now 
has to vanish on both sublattices for the zigzag termina- 
tion too. For the armchair ribbon, the superposition in 
Eq. (|6| consisted of plane waves around the K and 
points (as in the absence of spin-orbit coupling). In the 
zigzag case we need a superposition of two plane waves 
with different decay lengths in the transverse direction. 
These decay lengths, up to second order in e = E/t and 
r = tso/t, are given by A^ = 1/K.i = -^[log(^i)]"\ 
where 



2 Qi 



Qc 



(9a) 



Qc{Ql 



Here we have defined Qc = 2 cos(>/3gaja/2), Qi = 
2 sin(>/3gaja/2), and Q2 = 2 shi{-\/?>qxa) . The energy 
dispersion at the lower boundary is given by (see Ap- 
pendix) E = —6tsoSz sm{^/3qxa). The spin-up and spin- 
down bands thus cross at the wave vector QxCl = 7r/>/3 
with zero energy. The spin-up (down) electrons therefore 
propagate to the right along the lower (upper) boundary. 

At low energies and for small values of the spin-orbit 
coupling, the contribution ^2 in Eqs. ^ is clearly dom- 
inating. Furthermore, the spin-orbit coupling r only 
contributes a second order correction to (bo. We thus 
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FIG. 5: (Color online) Upper panel: Conductance vs Fermi 
energy for a QSH insulator/metal junction in an armchair 
ribbon with edge roughness of different depths d. The length 
of the disordered region is L = 100a. Every data point is 
obtained by averaging over 256 disorder configurations. Other 
parameters are as in Fig. [2] In the inset the same conductance 
is shown in the absence of spin orbit coupling. Lower panel: 
Ribbon geometry (finite section) for various roughness depths. 



find that the edge state width for zigzag ribbons is only 
weakly dependent on the exact value of the spin-orbit 
coupling. On the other hand, it is rather strongly de- 
pendent on the wavevector (through Qc), and thus 
on the energy. In Fig. |4] we have plotted the edge state 
widths A = for armchair and zigzag orientations as 
a function of the Fermi energy E. Except for energy val- 
ues close to the spin-orbit gap, the width of the zigzag 
edge states is considerably smaller than the armchair one. 
This is precisely what we found in Fig. [3j However, for 
energies close to Aso, the zigzag edge state width in- 
creases dramatically (as is also evident in Fig. 3|. For 
these energies, A also becomes dependent on the spin- 
orbit coupling strength, as shown in Fig. [4] with the red 
solid, dashed and dotted lines. 

To understand the perfect transmission found in zigzag 
ribbons, first consider a ribbon without spin-orbit cou- 
pling. In that case, one can perform the low energy ap- 
proximation around the two inequivalent ^ points in 
the same fashion as with armchair terminated ones. In 
the valley isotropic basis, the isospin is then simply given 
by = (1, 0) for a left moving particle and |r)^ = (0, 1) 
for a right moving one at the lower and upper boundaries. 
Edge states along zigzag boundaries thus belong to op- 
posite (orthogonal) valleys for opposite propagation di- 
rections. By turning on the spin-orbit coupling, the low 
energy dispersionless bands evolve into dispersive heli- 



FIG. 6: (Color online) The same as Fig. [sjbut for a zigzag 
terminated ribbon. 



cal bands crossing at the M point, as shown in Fig. [T] 
and their isospin is not well defined anymore. The lat- 
ter bands are no longer separated by a large momen- 
tum, and one might therefore think that the potential 
step, even if it is smooth, could induce scattering be- 
tween these bands. The key point now is that for a fixed 
spin, these two bands belong to opposite boundaries in 
the left region, so no large momentum can connect them. 
Scattering between these bands is only possible if the 
wavefunction is able to spread through the bulk of the 
ribbon. When this happens, to the right of the junc- 
tion, the Fermi energy of the quasiparticles lies above the 
bulk gap where the high energy modes in this region have 
again a well defined valley polarization. Hence, all modes 
injected from a left moving helical edge state will belong 
to the K valley and not to the K' one. The backscat- 
tering channel, that is derived from the valley, will 
therefore remain inaccessible. 



III. INFLUENCE OF EDGE ROUGHNESS 

For a clean ribbon, perfect transmission was due to the 
orthogonality between equal-spin left- and right-moving 
edge channels, combined with valley-preserving scatter- 
ing throughout all the ribbon. Disorder in general will 
break the orthogonality between topological edge chan- 
nels and a decrease in the transmission through a QSH in- 
sulator/metal junction is thus to be expected. However, 
even if the disorder conserves the orthogonality of the 
topological channels, e.g., when disorder is only present 
in the right (metallic) region of the junction, it will in 
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general induce inter- valley scattering which will connect 
counterpropagating egde channels through the bulk of 
the metallic region. Note that there are no topological 
edge states running along the junction, since this would 
require that both sides of the junction be bulk insulators 
with different values of the topological invariant. 

To analyze the influence of disorder in more detail we 
have considered rough edge boundaries, which break the 
orthogonality between topological edge states. We model 
rough edges by randomly removing carbon atoms up to 
a depth d (in units of a) from the ribbon's boundaries. 
In particular, we model a rough edge with a function of 
the form 

f[x) = Vh ^^-J^^-^'^^^X - ^n^ . (10) 

Here 7/5 is VK— (i/2 for the upper boundary and (i/2 for the 
lower one. The rough ribbon is obtained by removing all 
the carbon atoms in positions {x^y) for which y > fix) 
{y < fix)) at the upper (lower) boundary. The number 
of modes N in the Fourier series is chosen to be of the 
order of the number of atoms along the length of the 
ribbon. The amplitudes and phases (fn are random 
numbers between and 1 or — tt and tt respectively. Note 
that this function generates a boundary profile that varies 
strongly at the scale of the lattice spacing a and therefore 
is a source of strong inter- valley scattering. Examples of 
ribbons with different roughness depths are sketched in 
the lower panels of Figs. |5] and [6) where sections of width 
W = 54a (as used in the calculations) are shown. 

The conductance of a QSH insulator /metal junction in 
the presence of edge roughness is depicted in Figs. [5] and |6] 
for armchair and zigzag ribbons respectively. Results for 
intermediate ribbons are similar to one or the other (de- 
pending on the crystallographic orientation) and are not 
shown. All conductance points in these plots are ob- 
tained by averaging over 256 disorder realizations. The 
non-perfect transmission of the topological edge states 
through the junction is manifested by a reduction of 
the first conductance plateau. This reduction is however 
much less pronounced than the conductance decrease in 
a disordered n-n (or p-p) junction in the absence of spin- 
orbit coupling (and thus without topological edge states 
in the left region), as shown in the insets of Figs, [s] and |6] 

This increased protection against backscattering can 
be understood from plots of the charge density in the 
sample, as depicted in Fig. [7| In the presence of edge 
roughness, the spin-polarized edge state incoming from 
the upper boundary does not follow a straight line any- 
more (like in the flakes of Fig. |3| but its wavefuction 
accommodates to the roughened boundary. Therefore, 
unless d is so large that upper and lower topological edge 
state wavefunctions can overlap, there are still no avail- 
able intra-boundary backscattering channels. Only a sin- 
gle reflection channel is open, namely to cross to the topo- 
logical edge state at the lower boundary. This is radically 
different in a ribbon without intrinsic spin-orbit coupling. 



Armchair Zigzag 




FIG. 7: (Color online) Numerical results for the charge den- 
sity across the same potential step as in Fig. [3] but in the 
presence of edge disorder of depth d = 5a. The topological 
edge channel from the upper left boundary is transmitted into 
bulk normal channels to the right but can now be backscat- 
tered towards the lower left boundary. Only one type of spin 
current is considered (a similar contribution with the opposite 
spin, not shown, is incoming from the lower boundary). 

where multiple transmission and reflection channels are 
available throughout the bulk and disorder will have a 
bigger impact. The previous argument explains the small 
sensitivity of the conductance on roughness depth d for 
energies below the spin-orbit induced gap in Figs. |5] and 
|6] However, when the edge state width A becomes com- 
parable to the ribbon's effective width (given approxi- 
mately hy W — 2(i), direct inter-boundary backscatter- 
ing can take place in the left region due to the over- 
lap of counter-propagating (non-orthogonal) edge chan- 
nels, significantly reducing the transmission probability. 
This happens for zigzag ribbons with energies close to the 
band-gap edge, where the edge state width diverges (see 
Fig. [4]), which manifests itself as a conductance dip in 
Fig. |6[ For armchair ribbons, this dip does not appear as 
the edge state width is independent of energy and much 
smaller than the ribbon width for the parameters we have 
chosen. 

In Fig. [S] we plot the conductance through a QSH insu- 
lator/metal junction as a function of the ribbon's width 
W in the presence of edge roughness with a depth d = 5a. 
The disordered region has a length L = 100a, and the 
potential step parameters are Uq = 0.25t and S = 5a. 
Zigzag and armchair ribbons are considered, and ener- 
gies either close to the Dirac point or close to the spin- 
orbit induced gap are compared. For all parameters, we 
clearly see that the conductance increases with so 
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FIG. 8: (Color online) Conductance vs Fermi energy for a 
QSH insulator/metal n-n junction as a function of the rib- 
bon's width W. Edge disorder is present with roughness 
depth d = 5a and each conductance point is averaged over 
180 configurations. Other parameters are as in Figs. |5] and [6] 
Zigzag and armchair ribbons with energies either close to the 
Dirac point or to the topological band-gap edge are compared. 



that perfect transmission is recovered in the thermody- 
namic limit. This can be understood by realizing that, 
as the ribbon width is increased, the number of avail- 
able forward scattering channels grows while the avail- 
able backscattering channels remains one. Note that for 
the potential step parameters selected, all curves behave 
similarly except for the one of the zigzag terminated rib- 
bon for energies close to the gap. The slower conductance 
increase in that case is a result of the large overlap of the 
opposite boundary topological edge states, as discussed 
previously. Although the junction used for the calcula- 
tions in Fig. [S] is a n-n junction, we have found similar 
results for p-p and n-p (p-n) junctions and conclude that 
conductance recovery with increasing ribbon width is a 
general trend. 



IV. CONCLUSIONS 

Transport through QSH insulator /metal junctions has 
been analyzed in graphene ribbons in the presence of an 
electrostatic potential step. We have found that for clean 
ballistic ribbons (with arbitrary crystallographic orienta- 
tion), the topological edge states are transmitted per- 
fectly through such a junction, independent of the rib- 
bon's width, the parameters of the potential step and 
the value of the intrinsic spin-orbit coupling, as long as 
inter- valley scattering is not present. This has been ex- 
plained by proving the isospin orthogonality between for- 
ward and backward edge channels analytically, a mech- 
anism that is analog to the Klein tunneling. We have 
also found that the edge state width only depends on the 
spin-orbit coupling strength for armchair ribbons but is 
in addition strongly energy dependent for zigzag ribbons. 
In the presence of edge roughness, the orthogonality, and 
thereby the perfect transmission, is destroyed and the 



upper (lower) edge channel can be reflected to the lower 
(upper) one through the metallic region where the parti- 
cle's wave function spans the whole ribbon's width. How- 
ever, even in the presence of strong disorder, reflection 
from a QSH insulator/metal junction is small as com- 
pared to a similar ribbon in the absence of intrinsic spin- 
orbit coupling. This protection can be understood as a 
spin-orbit proximity effect^^ since to the left of the junc- 
tion the edge states are still topologically protected (as 
long as time-reversal symmetry breaking disorder is ab- 
sent). Finally, we found that perfect transmission is grad- 
ually restored in the presence of disorder as the ribbon's 
width increases, since the number of available forward- 
scattering channels grows. This recovery is dependent 
on the ribbon's crystallographic orientation, and was ex- 
plained by analytical considerations as well as by numer- 
ical simulations of the charge density across the junction. 
In the present work, the QSH insulator/metal junction 
was realized by means of a constant spin-orbit coupling 
and an electrostatic potential step. An analogous junc- 
tion could be achieved in the presence of a spin-orbit 
potential step if the intrinsic spin-orbit coupling could 
be externally tuned^^ across the ribbon. We analyzed 
this scenario numerically and found similar results for 
metallic ribbons. 



Acknowledgments 

EP acknowledges funding from the MICINN- Spain via 
grant FIS2009-08744. 



Appendix A: Topological Edge States 

In this Appendix we derive the edge state energies and 
wavefunctions for armchair and zigzag graphene ribbons 
in the presence of the intrinsic spin-orbit interaction. 
Moreover, we prove the orthogonality between same spin 
left and right moving edge states at the opposite bound- 
aries of the ribbon and obtain the edge state width as a 
function of spin-orbit coupling and energy. 

We consider the lattice and the reference system shown 
in Fig. [9] In terms of the lattice vectors bi (linking 
nearest-neighbor carbon atoms) and (linking next- 
nearest neighbors), the tight-binding Hamiltonian can be 
written diS H = Hq + Hso^ where 



ri+bn 



h.c, 



(Al) 



and 



n=l 



it 



ffso = itso^^(-ir+'s. 

~bi,s. bf,+s„,s. ] ■ (A2) 
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FIG. 9: Geometry of the graphene lattice considered in 
the derivation of the edge state wavefunctions. Here 61,2 = 
a(±V3/2,0), h = a(0,3/2), ai,2 = a(=bV3/2, 3/2), a6,3 = 
a(div^, 0) and as, 4 = a(div^/2, — 3/2), where a is the inter- 
atomic distance. 



In the previous two equations, i runs over equal sublattice 
sites, and are creation operators on sublattice A and 
B respectively, and 5^ = ±1 is the z-component of the 
electron spin. 

Fourier transforming these Hamiltonians into momen- 
tum space, one gets: 



Ho = Yl 4,5. ^^^^^ + 



(A3) 



where q is the wave vector in the graphene plane mea- 
sured from the center of the first Brillouin zone and 



(A4) 

with a the interatomic distance. On the other hand 



with 



X{s.,q) = is Jso Yl (-ir+'e-'«-" 



(A6) 



n=l,...,6 



3 

2sin(^a(7^)cos(-a(7^) - sm{V3aqx) 



The total Hamiltonian can thus be written as 



with eigenvalues 



E{q) = ±Vm\' + x{s,,q)' 



(A7) 
(A8) 



and eigenvectors 



(X 



The total wave function is then 



1 



(A9) 



^.^,,Kx,,)=e--e--(j-j^;;J^^ (AlO) 



Armchair ribbon 



With the geometry depicted in Fig. [9j armchair edges 
extend along the ^-direction. To obtain an edge state 
localized in the x-direction, we replace Qx i/^ in 



Eq. (AlO). The boundary condition amounts to setting 
the wavefunction to zero at x = 0, both on A and B 
sublattice sites. To obtain a nontrivial solution, we have 
to consider a superposition of two wavefunctions. In the 
armchair case, one can consider a combination of waves 
around the K and K' points of the Brillouin zone of 
graphene, because the low-energy limit corresponds to 
having wavevectors in the vicinity of these points. Our 
ansatz is thus 



ao e 



(All) 

where Qy is measured from the K'^ ^ point (with K'^ ^ = 
±^^^x) and ao^Po are the relative amplitudes of the 
wavefunction in each valley. 



After expanding the quantities ^ and x Eqs. (A4) 
and (A6) around the K^''' point, one arrives at ^ = 
hvF{^qx + iQy) and ^ = i'5;2Aso7 where the upper 
(lower) sign corresponds to t he K {K') point, and Aso = 
3A/3tso- The spinors in Eq. ( All ) are then given by (c.f. 
Eq. 



(A12) 



A nontrivial wavefunction (All) fulfilling the bound- 



ary condition ^^^(x = 0, y) = ^ ^ ^ can only be found 

(A13) 



when 



i^A V^A 



0. 



This leads to the following expression for K: = 
—Sz ^y^^ . Inserting this into the dispersion relation E = 



± A'^Q -\- {hv f)'^ {q^ — K.'^) ^ one obtains E = ±.hvFqy 
The result for n can then be written as: 



= Sz sgn(^^) sgn(^). 

flVF 



(A14) 
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At the left boundary, n should be positive to get a decay- 
ing solution, so, e.g., electrons {E > 0) propagate either 
with spin- up downwards {qy < 0), or with spin-down 
upwards Qy > 0. The pseudospin associated with these 
electrons is 



^9 
^9 



1 

71 



1 



(A15) 



At the right boundary one needs n <{) and the propaga- 
tion directions are reversed. The pseudospin is then 




1 

7i 



5^1 
1 



(A16) 



In terms of the Pauli matrices for the pseudospin de- 
gree of freedom, a, and for the isospin one, r, the Hamil- 
tonian (A7) around the K and K' valleys can be ex- 
pressed: 

H = hvriqxcrx ^ Tz - Qyay ^ To) ^ SzAsoa^ ^ Tz, (A17) 

acting on the Dirac 4-component spinor tj; = 
(V^A, V^B, '0^5 '^b)^- order to prove the orthogonality 
between equal-spin incoming and outgoing edge channels, 
it is convenient to express the previous Hamiltonian in 
the "valley isotropic" basis,'^ where H is isotropic in the 
valley degree of freedom. The Hamiltonian in the val- 
ley isotropic basis is related to H through the unitary 
transformation = U^HU, where 



= [hvriqxcrx - QyCTy) + s^Asocr^] <S) tq (A18) 



and 



U = 2(^0 +r,) 



1 



) (Jo 



(to - Tz) 



Jy. 



(A19) 



In a similar way, the pseudospins are transformed like 
^iso _ — {i^A,i^B,i^B, ~^'a)^ • The total wave func- 
tion in this basis is 



ae 



(A20) 

where a and {3 represent the valley polaritazion or isospin 
of the wave function. 



2. Zigzag ribbon 

Now we consider a semi-infinite zigzag-terminated 
graphene sheet to find expressions for the topological 
edge states. In the absence of spin-orbit coupling, it is 
known that the boundary condition amounts to a van- 
ishing wavefunction on one of the sublattices. However, 
as spin-orbit coupling induces next-nearest neighbor hop- 
ping, the total wave function in our case has to vanish 
in both the A and the B sublattice sites at the edge. 



This condition can never be met by a single wavefunc- 
tion of the form (AlO), so we will look at a superposition 
of two of them. Furthermore, as we are looking for an 
edge state solution the wavefunction has to decay from 
the edge into the bulk. Considering the lower boundary 
(the semi-infinite graphene sheet extends towards ^ > 0, 
see Fig. [9|, we then substitute qy m and make the 
following ansatz: 



a e 



+ ^ e- 



In the following, the notation 



i^Aisz.qx.i^i) 
i^B{sz,qx,i^i) 

'^A{Sz,qx,l^2) 
'^B{Sz,qx-,l^2) 



.(A21) 



Qc 


2cos( — g^a). 


(A22a) 


Qi 


= 2sm(^(7^a), 


(A22b) 


Q2 


= 2 sin(A/3g'a^a), 


(A22c) 


(t>i 




(A22d) 



will prove to be convenient. The terms ^ and x 
Eqs. (A4) and (A6) can then be written as 



X{Sz,qxAi) = -tsOSz [Q2 - Ql{(l)i 

and the energy dispersion (|A8[) becomes 



{Qc 



{Qc 



(A23a) 
-^)](A23b) 



H)^r' (g2-Qi(</>z + 

(A24) 

Here we defined e = E/t, and r = tso/t. Multiplying 
this equation with 0f gives a fourth order equation in 
from which one can in principle calculate the allowed 
values of as a function of energy. Solving this equation 
exactly proves to be too difficult. However, as we are 
interested in small values of the energy e and spin-orbit 
coupling strength r, we can make an expansion to second 
order in these parameters. From the four solutions to 



'Eq. (A24), only two fulfill the condition 1 < 1 for a 
decaying wavefunction around the M-point {qxa = 7r/\/3 
corresponding to zero energy e = 0). They are 



>^Ql 
Qc' 



Qc 



1 



> (Qc(Q2 + QcQi) + Qi)^ 
Qc{Q'c - 1) 



(A25a) 



(A25b) 



The quantities still depend explicitly on both en- 
ergy and wavevector qx. In order to find the dispersion 
relation e{qx)^ we have to apply the boundary condition 





ix,y = 0) 



to the wavefunction defined in 
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and /3 is then obtained when 



Eq. (A21). A non-trivial solution for the coefficients a 

= 0. (A26) 



With Eq. ( |A9[ ) for the spinors, this gives 

^{Qx^ h)[E{sz,qx) - X{sz^ Qx^ ^i)]. 



(A27) 



Using Eqs. ( |A23|) and inserting the allowed values for the 
(pi from Eq. ( A25 ), one obtains the dispersion relation (up 



to ffist order in r): 



^{sz.Qx) = -^SzrQ2 = -6sz^^ sin{V3qxa). (A28) 

The zero of energy lies at g^^a = tt/VS (at g^^a = there 
is no edge state solution with \<pi\ < 1). Quasiparti- 
cles with spin-up (spin-down) thus move with positive 
(negative) velocity at the lower boundary. At the 

upper boundary, the role of spin-up and spin-down is in- 
terchanged (see below). 

The in Eq. (A25) can assume negative values. With 
the definition 0^ = e~^^'^% this means that the solutions 
for can be complex and the wavefunction oscillates in 
the ^/-direction with a decaying envelope. In particular, 
we find that the wavefunction changes sign when propa- 
gating from one zigzag dimer line to the next. The edge 



state width has thus to be defined as = l/Re{hZi). 



With Eq. (A22d), one obtains 



3a 



21og(|(/>.|)* 



(A29) 



Until now, we only considered the lower boundary. For 
the upper boundary, the semi-infinite graphene sheet ex- 
tends now to negative ^/-values. A localized state is then 
found by substituting qy 



-Qy 



This will interchange the off-diagonal terms ^ and ^'^ m 



iti in Eq. (|ATO). 

— 



the Hamiltonian ( A7). If we further substitute Sz —Sz-, 
then the diagonal terms are interchanged too. The 
Hamiltonian in Eq. (A7) for the upper boundary then 



looks exactly as the one for the lower boundary, but with 
the A- and ^-si tes interchanged. Using Eq. ( A9) and the 
condition (A27), it is possible to see that the pseudospin 



for the lower boundary is 



W{sz,qx,(t>i)) = 



SzTQi 



(A30) 



to ffist order in r, both for z = 1,2. With the substitution 
A ^ 5 and Sz —s^ mentioned above, the pseudospin 
for the upper boundary is 



W^Sz^qxAi)) 



1 



(A31) 
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